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1  Introduction 


PDRC  Report  86-09  discussed  an  averaging  algorithm  to  solve  the  aggregate- 
disaggregate  reformulation  of  MRMATE.  The  procedure  was  presented  as 
an  alternative  to  Benders’  decomposition  for  solving  the  LIFTCAP/MRMATE 
i  linking  problems  in  MODES  in  PDRC  Report  87-06.  This  report  presents 

two  approaches  for  solving  MODES. 

i  The  first  approach  formulates  MODES  as  a  hierarchy  of  decisions  to  be 

generated.  The  hierarchy  consists  of  an  aggregated  version  of  LIFTCAP 
\  that  generates  allocation  of  assets  to  aggregate  channels.  An  aggregate 

MRMATE  model  allocates  MRs  to  these  aggregated  channels  generated 
f  by  aggregate  LIFTCAP.  An  aggregate  LIFTCAP  model  and  an  aggregate 

MRMATE  model  constitute  the  aggregate  MODES  model. 

'  A  solution  to  the  aggregate  MODES  model  is  then  disaggregated  by  dis- 

r  aggregate  LIFTCAP  and  disaggregate  MRMATE  model.  The  disaggregate 

problems  generate  solutions  at  the  required  level  of  detail  consistent  with 
aggregate  MODES  solutions.  Consistency  implies  that  detailed  channel  ca¬ 
pabilities  and  MR  allocations  added  across  an  aggregate  bundle  is  bounded 
by  the  aggregate  channel  capabilities  and  MR  allocations  generated  by  the 
aggregate  MODES  model. 

The  second  approach  considers  the  MODES  problem  from  a  different 
perspective.  In  some  situations,  the  dominant  objective  is  to  move  all  MRs 
to  their  destination  within  their  desired  time  window  using  the  minimum 
number  of  assets.  The  driving  problem  in  this  version  of  modes  is  the 
MRMATE  problem.  MRMATE  generates  allocations  of  MRs  to  their  best 
’  (minimum  penalty)  channels.  This  defines  channel  capabilities. 

Given  channel  capabilities,  the  corresponding  assets  are  required  to  re¬ 
alize  these  channel  capabilities.  This  yields  a  model  MIN  ASSET  which 
solves  the  problem  of  generating  the  minimum  number  of  assets  required. 
MINASSET  also  generates  a  cost  of  using  a  channel  for  the  MRs.  Interac¬ 
tion  between  MRMATE  and  MINASSET  determines  an  allocation  of  MRs 
to  channels  that  minimizes  assets  required  to  affect  the  delivery  pattern 
provided  as  input  to  the  MODES  model. 
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2  A  Hierarchical  approach  to  MODES 

A  desirable  feature  of  a  hierarchical  version  of  MODES  is  separation  of 
the  decision  making  process.  The  effect  is  that  solutions  generated  by 
the  procedure  are  consistent  with  desirable  aggregate  MODES  solutions. 
Since  the  aggregate  MODES  model  is  of  a  smaller  size,  many  iterations 
of  aggregate  MODES  would  be  run  until  acceptable  aggregate  decisions 
are  generated,  then  these  decisions  would  be  disaggregated  to  generate  a 
MODES  solution.  Furthermore,  aggregate  solutions  enable  decomposition 
of  the  disaggregation  process  by  region.  This  allows  regional  problems  to 
be  solved  separately. 

2.1  Aggregation  of  ports 

The  source  region  of  POEs  and  the  destination  region  of  PODs  are  divided 
into  regions.  All  channels  between  a  pair  of  source-  destination  regions  have 
the  same  cycle  time.  This  implies  that  assets  operating  between  POE-POD 
pairs  in  a  pair  of  regions  have  the  same  travel  time. 

The  aggregate  LIFTCAP  model  consists  of  a  network  flow  problem  for 
channels  between  aggregate  PODs  and  POEs  and  asset  side  constraints 
which  limit  total  asset  allocation  across  aggregate  channels.  Given  an  ag¬ 
gregate  LIFTCAP  solution,  the  disaggregate  problems  create  one  disaggre¬ 
gate  LIFTCAP  problem  for  each  aggregate  POE-POD  asset  triple. 

LIFTCAP  POEs  are  divided  into  L  sets  Pi,...,  Pi  and  PODs  are  di¬ 
vided  into  M  sets  S\,. . .  ,Sm-  The  aggregate  LIFTCAP  model  is: 

Ea  E,  sija  <Zk(p,Ck  for  i  =  1,2 ,...L 
Ea  E<  Sijo  <  EieS;  Dk  for  j  =  1,2, ...  Af 
Ei  Ej  ®«ja  5;  Aa  for  a  =  1,2,. . .  B 

Sija  ^  0 

Given  an  aggregate  LIFTCAP  solution,  the  disaggregate  LIFTCAP 
problems  are  created  as  follows.  For  each  of  the  LxMxB  sets  of  aggre¬ 
gated  channels,  the  disaggregation  problem  Pi  —  Sj  -  a  generates  detailed 
channel  capabilities.  Disaggregation  problem  P,  —  S;  -  a  is: 
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£  k  t  S3  clka  <  C{  for  1  e  P, 

LktP,  ckma  <  Dj  for  m  e  S , 

£l  £m  clma  —  ^ija 

The  aggregate  and  disaggregate  LIFTCAP  models  have  the  matrix 
structure  in  Figure  1. 


r- — 

AGGREGATE 

LIFTCAP 

■fa  4].  4:.. 

h  n  r 

=3  ..  <  SiJa 

_ ll JO 

Pj-Sj-a 

— 

DISAGGREGATE 

LIFTCAP 


Figure  1:  Aggregate  and  Disaggregate  LIFTCAP  structure 


2.2  MRMATE  Aggregation 

PDRC  Report  86-04  discussed  different  functional  aggregation  strategies 
for  solution  of  MRMATE.  The  objective  here  is  to  create  an  aggregate 
MRMATE  model  that  allocates  MRs  to  aggregate  channels.  The  aggregate 
channel  capabilities  are  disaggregated  by  disaggregation  problems  which 
generate  detailed  MR  allocations  to  detailed  channels. 
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Channels  are  aggregated  into  sets  of  LxMxB  aggregate  channels  with 
capabilities  provided  by  aggregate  LIFTCAP  (  i.e.  s„a) .  Thus  the  aggre¬ 
gate  MRMATE  problem  is  expressed  as 

Zk<F(i.j,a)  yrk  =  Mr  for  r  =  1,2,...  T 
Er  Vrk  <  Sija  f  OT  k  =  1,2,...  ( LlMxB ) 

Vrk  >  0 


MRMATE  Disaggregation  problems  assign  MRs  to  detailed  channels. 
This  results  in  a  disaggregation  problem  for  each  set  of  aggregated  channels, 
yielding  in  LxMxB  disaggregation  problems. 

Disaggregation  problem  ija  becomes 

min  El  Em  Plm  Xlm 

s.t.  m  €  F(i,j,a)  xlm  =  ylk  for  1  =  1,2, ...  T 
Ei  Z(m  =  Ck  for  k  €  F(i,j,a) 

Xlm  >  0 

The  aggregate  and  disaggregate  MRMATE  models  have  the  structure 
in  Figure  2. 

2.3  Aggregated  LIFTCAP  and  MRMATE  Models 

Consider  the  aggregate  LIFTCAP  and  MRMATE  models  in  conjunction. 
It  has  the  structure  in  Figure  3.  The  aggregate  LIFTCAP  model  generates 
channel  capabilities  required  at  the  level  of  the  aggregate  MRMATE  model. 
Rearranging  the  blocks  in  Figure  3  a  structure  as  in  Figure  4  is  produced. 

Figure  4  presents  a  hierarchical  model  of  MODES.  Aggregate  LIFTCAP 
and  aggregate  MRMATE  models  represent  the  aggregate  MODES  models. 
The  aggregate  MODES  model  generates  aggregate  channel  capabilities  and 
an  allocation  of  MRs  to  aggregate  channels.  Aggregate  channels  represent 
an  allocation  of  MRs  to  regions,  i.e.  POE-POD  regions.  The  aggregate 
MODES  model  has  the  same  structure  as  the  MODES  model  but  is  of  a 
much  smaller  size. 

Smaller  size  of  the  aggregate  MODES  model  implies  that  various  asset 
allocation  scenarios  can  be  analyzed  to  determine  a  suitable  aggregate  asset 
allocation  very  quickly. 
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AGGREGATE  MODES  MODEL 


Figure  4:  Hierarchical  MODES  Model 


'1 


The  disaggregate  MODES  model  attempts  to  generate  asset  allocations 
and  channel  capabilities  consistent  with  the  aggregate  solution.  This  im¬ 
plies  that  each  POE-POD  region  can  be  treated  independently  in  the  dis¬ 
aggregate  MODES  model  to  generate  MR  allocations  at  each  regional  level 

consistent  with  aggregate  asset  allocations  and  MR  assignments.  | 

The  disaggregate  MODES  model  generates  detailed  channel  capabilities  i 

consistent  with  aggregate  asset  allocations  and  detailed  MR  allocations  to  j 

these  channels.  Given  an  aggregate  solution,  from  in  Figure  4  it  can  be  ' 

seen  that  the  disaggregate  problems  separate  into  LxMxB  disaggregate  I 

problems,  one  for  each  POE-POD  region  and  asset  type.  The  structure  of  \ 

the  MODES  disaggregate  problems,  for  a  given  aggregate  solution,  is  the  1 

same  as  the  MODES  structure  itself. 

The  disaggregate  problems  are  again  of  much  smaller  size  than  the  orig-  ' 

inal  MODES  model.  Furthermore,  the  disaggregate  problems  are  restricted  | 

to  each  POD-POE  region  and  asset  type.  This  implies  that  the  disaggre-  I 

gate  problems  can  also  be  used  to  generate  detailed  problem  information  j 

in  terms  of  port  capacity  and  channel  capability. 

2.4  A  solution  procedure  using  averaging 

The  problem  structure  in  Figure  4  could  be  exploited  the  aggregate  and 
disaggregate  MODES  problems  could  be  decomposed  in  such  a  manner  as  to 
maintain  structure  at  each  iteration  of  an  iterative  solution  process.  PDRC 
Report  86-09  discusses  use  of  an  averaging  procedure  to  solve  re-formulated 
aggregate  disaggregate  MRMATE.  The  algorithm  is  asymptotically  optimal 
and  provides  the  desirable  features  outlined  above. 

The  algorithm  would  work  as  follows: 

1.  Initialize  the  solution  process  by  generating  an  aggregate  MODES 
solution  that  is  disaggregated  to  generate  duals  for  the  disaggregation 
problems. 

2.  Use  the  disaggregate  problem  duals,  associated  with  each  of  the  ag¬ 
gregate  channels  and  aggregate  POE-POD  arcs,  to  generate  an  ag¬ 
gregate  MODES  solution  optimal  with  respect  to  these  costs.  (The 
aggregate  problem  objective  function  value  generates  a  lower  bound 
to  the  overall  MODES  problem.) 
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3.  Average  the  aggregate  solutions  generated  up  to  and  including  this 
iteration  and  disaggregate  the  averaged  aggregate  solution  generated. 

4.  The  aggregate  solution,  when  disaggregated,  generates  disaggregate 
duals  which  are  again  averaged  across  all  disaggregate  duals  generated 
through  this  iteration.  (The  disaggregate  problems  generate  an  upper 
bound  to  the  overall  MODES  problem.) 

5.  Repeat  steps  2  through  4  until  the  upper  and  lower  bounds  differ  by 
less  than  some  prespecified  value. 

3  An  alternative  approach  to  MODES 

This  approach,  called  ALLOCATE,  considers  relaxation  of  asset  constraints 
and  port  capability  constraints  at  each  iteration  of  the  model.  Successive 
iterations  improve  asset  and  port  capability  feasibility  of  MR  allocations 
generated.  At  every  iteration,  delivery  patterns  are  generated  which  fall 
within  specified  time  windows.  If  the  assets  available  are  less  than  required 
to  affect  the  specified  delivery  pattern,  then  the  procedure  generates  the 
minimum  assets  required  to  maintain  this  delivery  pattern. 

3.1  Allocation  of  MRs  to  channels 

The  master  problem  is  the  allocation  of  MRs  to  channels.  This  is  a  trans¬ 
portation  problem  with  channel  capabilities,  i.e.  supplies  unspecified  and 
costs  depending  on  channel  capability  usage.  The  objective  is  to  generate 
an  allocation  of  MRs  to  channels,  which  specifies  the  channel  capabilities, 
that  minimizes  MR  to  channel  allocation  penalties  as  well  as  generates  the 
minimum  cost  channel  allocation. 

min  Si  Sy  Pij  %ij  4"  Sy  Cy 
s.t.  x ij  =  m,  for  i  =  1,2 ,...,M 

Si  *«y  -  Cj  =  0  for  j  =  i,2,...,JV 
Xij  >  0 

The  optimal  solution  to  this  problem  would  result  in  each  MR  choosing 
a  single  channel.  This  is  a  desirable  feature  if  MR  split  across  channels  is 
not  desirable. 


11 


3.2  Determining  Asset  requirements 

Given  a  set  of  channel  capabilities  generated  by  the  master  problem  of 
ALLOCATE,  the  minimum  number  of  assets  required  to  realize  that  set 
of  capabilities  must  be  determined.  This  problem  is  called  MINASSET. 
Furthermore,  if  a  set  of  costs  associated  with  each  of  the  asset  types  along 
with  capability  constraints  is  given,  the  subproblem  determines  a  minimum 
set  of  assets  that  realize  the  capabilities  specified. 

The  minimal  number  of  assets  is  determined  by  creating  a  minimum 
cost  network  flow  model.  Nodes  in  the  model  represent  channels.  An  arc  is 
included  from  node  (channel)  i  to  node  (channel)  j  if  an  asset  can  complete 
the  task  of  moving  MRs  along  channel  i  and  travel  to  the  POE  of  channel 
j  and  get  to  its  destination  before  the  time  specified  by  channel  j. 

Each  asset  type  generates  a  different  minimal  cost  network.  Asset  ca¬ 
pacities  provide  the  number  of  units  of  each  asset  required  to  realize  channel 
capabilities. 


3.3  Generating  channel  capability  costs 

When  the  subproblem  is  solved,  a  set  of  asset  requirements  is  generated 
along  with  a  dual  variable  for  each  of  the  nodes  (channels).  Duals  associated 
with  nodes  reflect  the  cost  of  using  that  channel  for  any  MR.  These  duals 
are  used  as  costs  on  channels  in  the  master  problem  of  ALLOCATE.  New 
costs  are  used  to  generate  a  revised  MR  allocation  with  respect  to  assets 
required  and  this  completes  an  iteration. 

3.4  A  solution  procedure  using  averaging 

The  solution  procedure  is  based  on  the  averaging  approach  described  in 
PDRC  Report  86-09.  It  has  the  property  that  problem  structure  is  main¬ 
tained  for  both  the  master  and  subproblems  at  each  iteration.  It  converges 
asymptotically  to  the  global  optimal  solution.  The  algorithm  works  as  fol¬ 
lows: 

1.  Solve  MRMATE  and  generate  MR  allocations  to  channels.  This  de¬ 
termines  a  set  of  node  capability  constraints  for  the  subproblem. 
Channel  costs  in  the  first  iteration  are  set  to  some  reasonable  es¬ 
timates. 
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2.  For  the  channel  capabilities  generated,  solve  the  subproblem  to  gen¬ 
erate  a  minimal  number  of  assets  required,  together  with  duals  for  all 
channels. 

3.  Solve  the  master  problem  with  these  duals  as  costs  on  channels  and 
generate  a  new  MR  allocation.  (This  generates  lower  bounds  to  the 
overall  problem.) 

4.  Average  across  all  MR  allocations  generated  thus  far  and  use  that 
average  solution  for  the  subproblem. 

5.  Solve  the  subproblem  with  this  new  MR  allocation  to  generate  chan¬ 
nel  duals.  (This  problem  generates  upper  bounds  on  the  objective 
function  value.) 

6.  Average  across  all  channel  duals  and  use  this  information  in  the  mas¬ 
ter  problem.  Stop  if  the  difference  between  upper  and  lower  bounds 
is  less  than  some  prespecified  value. 

4  Conclusions 

The  two  approaches  presented  in  this  report  present  alternatives  to  the 
current  implementation  of  MODES. 

The  hierarchical  MODES  model  attempts  to  reduce  problem  size  while 
maintaining  structure.  The  strategy  is  to  separate  out  MODES  problems 
by  POE-POD  regions  and  enable  independent  analysis.  It  also  enables 
various  alternative  asset  allocation  scenarios  to  be  tested  quickly  at  the 
aggregate  level  before  analyzing  the  problem  in  detail.  Finally,  since  the 
averaging  procedure  maintains  problem  structure  at  each  iteration,  efficient 
codes  can  be  used  to  solve  the  problem  at  each  iteration. 

ALLOCATE  represents  an  approach  in  which  the  emphasis  is  on  gener¬ 
ating  delivery  patterns  at  the  cost  of  the  availability  constraints.  Successive 
iterations  generate  delivery  patterns  which  use  a  lesser  number  of  assets. 
Stopping  at  any  iteration  would  provide  the  desired  delivery  pattern  but 
would  use  more  assets  than  available.  This  model  is  suitable  in  situations 
wherein  the  number  of  assets  must  also  to  be  determined  as  a  parameter 
in  the  model.  (Often,  the  asset  constraint  is  not  a  hard  constraint.)  Thus 
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if  the  procedure  is  stopped  at  a  solution  wherein  the  assets  required  do  not 
differ  considerably  from  the  number  of  assets  available,  then  a  satisfactory 
solution  has  been  obtained. 

The  disadvantage  of  ALLOCATE,  as  described,  is  that  it  ignores  port 
capability  constraints.  Procedures  are  required  which  would  incorporate 
this  information  into  the  ALLOCATE  model  while  maintaining  exploitable 
structure.  One  alternative  is  to  assume  that  port  capability  is  specified  by 
asset  type.  If  so,  then  this  information  can  be  included  in  the  form  of  a 
dummy  node  in  LIFTCAP. 


14 


